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ABSTRACT 

A method of predicting local values of heat transfer 
coefficients in round tubes for forced convection boiling 
of water with net vapor generation is proposed. The total 
heat transfer is postulated to be made up of a forced 
convection component and a nucleate boiling component 
where it exists. 

The forced convection compnnent was correlated by a 
Traviss (17) theoretical analysis of heat transfer across 
a thin annular film, modified to fit non-boiling data. 

The heat flux required to initiate nucleate boiling 
was predicted by the equilibrium of a hemispherical vapor 
bubble in a linear temperature gradient near the wall. 
The largest equivalent surface cavity radius capable of 
nucleation was suggested to be of order nos feet. 


Nucleate boiling correlations of Rohsenow (9), Mikic 





(22), and Thom (23) were examined to account for the boiling 
component. Superposition was accomplished by forcing the 
boiling component to be zero at the onset of nucleate 
boiling. 

The method was tested against eight sets of water 
data in vertical up and down flow. The Chen (8) correl- 
ation for convective boiling was also tested as a standard. 
The modified Traviss forced convection/Mikic nucleate 
boiling had the lowest average percent deviation between 
predicted and experimental values of wall superheat over 


all data points of +15.4%. 
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NOMENCLATURE 


Parameter used in incipient boiling analysis, ft?R 
Boiling number, q/A, 
hr, 


Constant in Mikic pool boiling analysis 


Specific heat capacity, BTU/lbm-°F 
Constant in Rohsenow nucleate boiling correlation 


Tube inside diameter, ft. 


Chen Reynolds number factor, (Re/Re, )9*9 


Traviss two-phase forced convection parameter, 
NuF 
2 


0.9 

1 Pri 

Traviss velocity profile parameter, Equation (2.4) 
Dengler boiling correction factor, Equation (1.2) 
Mass flow velocity based on total mass flow rate, 
lbm/hr-ft^ 

Acceleration of gravity, 4.173X10° ft/hr* 


8 ibm-ft/lbf- r^ 


Gravitational constant, 4.173X10 


Two-phase heat transfer coefficient, BTU/hr-ft~-°F 


Latent heat of vaporization, BTU/lbm 


Heat transfer coefficient obtained from Dittus- 
Boelter equation assuming total flow is all liquid, 


BTU/hr-ft^-?p 





(H-TS) 


Nu 


Pr 


3 u 


= < 


p$ 


tt 


Psat 


Gibbs free energy, BTU 
Thermal conductivity, BTU/hr-ft-?F 
Constant in Mikic pool boiling correlation 


Slope of fully developed boiling curve on log-log 
coordinates 


Nusselt number, hD 


Pressure, psia 


Prandtl number, “pA 
k 


Heat flux, BTU/hr-ft^-^F 


Radius of bubble cavity or bubble, ft 
Gas constant, ft-lbf/lbm-°R 


Chen two-phase effective Reynolds number 


Reynolds number for liquid fraction, GD(1-x) 
Mı 


Chen nucleate boiling suppression factor 
Temperature, °F or CR in difference equations, 
°R otherwise 


Specific volume, ft2/1bm 
Constant in Thom subcooled boiling correlation 


Vapor mass fraction, "quality" 


0.5 0.1 
Martinelli parameter, ES E 
x n 


My 
Perpindicular distance away from heated surface, ft 
Difference in vapor pressure corresponding to 


2 
A T ato lbf/ft 








AT 


A 
Tsat 


Chen effective superheat with flow, OF 


ο 
Wall superheat, ul); F 


Absolute viscosity, lbm/hr-ft 
Density lbm/ft? 
Surface tension, lbf/ft 


Parameter in Mikic pool boiling correlation 


Subscripts 


B 
crit 


data 


FC 


fg 
ib 


max 
pred 
sat 
ΠΡ 


try 


Fully developed nucleate boiling 
First cavity to nucleate 
Experimental value 

Effective value with flow 

Forced convection without boiling 
Associated with a change of phase 
Value at incipient nucleation point 
Value for liquid 

Largest cavity potentially active 
Predicted value 

Value at saturation conditions 
Two-phase 

Iterative trial value 

Value for vapor 


Evaluated at wall conditions 








Two-phase 
Iterative trial value 
Value for saturated vapor 


Evaluated at the heated wall 
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Chapter 1 


INTRODUCTION 


1.1 General 

There are countless applications where the transfer 
of heat from a surface to a boiling fluid is of great 
importance. Over the past 30 years many investigators 
have obtained experimental data, but attempts at describing 
the mechanisms and predicting heat transfer performance 
have usually met with limited success. The mass of cor- 
relations have in general had little applicability to 
systems much different from the one used to generate 
the data. 

It has become &lmost trite to cite the complexities 
of forced convection boiling phenomena, but since these are 
at the heart of the difficulty in predicting data, some 
of the considerations are reviewed in the following 
paragraphs. 

In single phase heat transfer, the heat flux generally 
varies linearly with the temperature difference, and is 
well predicted by any of several well known equations 
(Dittus-Boelter, Sieder-Tate, etc.). In nucleate pool 
boiling, however, the heat flux typically varies as the 
cube of the temperature difference, as shown in Figure 1. 


For non-boiling convection, the hydrodynamic field may be 
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FIGURE 1 Typical Curve for Saturated Pool Boiling 
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treated independently of the temperature field; but this is 
not the case for boiling phenomena (5). 
The case of convective boiling with net generation 
of vapor is shown in Figure 2. A vertical tube is heated 
uniformly over its length and is fed with subcooled liquid 
at a rate such that the liquid is completely evaporated. 
Experimental evidence indicates that several modes of 
heat transfer will occur as vaporization proceeds: 
1. Single phase forced convection to the liquid. 
2. Subcooled boiling. 
3. Saturated nucleate boiling. 
4, Two-phase forced convection (evaporation 
at the liquid film-vapor core interface). 
5. Combinations of modes 3 and L. 
6. Transition to a dry wall (liquid deficient). 
7. Dry wall (single phase forced convection 
to the vapor). 
In Figure 2 the fluid bulk temperature increases 
until saturation conditions are reached, and then gradually 
decreases with decreasing saturation pressure. The wall 
temperature profile initially increases parallel to the 
fluid temperature profile. At some point along the tube 
the conditions adjacent to the wall are such that bubbles 
of vapor can occur at nucleation sites. This mechanism is 
known as subcooled boiling. Because the liquid is subcooled 


(i.e. the bulk mean temperature is less than the saturation 
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temperature), any vapor bubbles that detach from the wall 
are condensed in the colder core liquid. As subcooled 
boiling begins, the wall temperature tends toward a constant 
value in accordance with the experimental observations 
that the heat flux in this region is a strong function 
of the wall superheat, Tw Tsat’ alone. 

The transition from subcooled to saturated nucleate 
boiling occurs when the thermodynamic mixed mean enthalpy 
equals the saturation enthalpy (x=0.). At this point, 
there is still subcooled liquid in the core, which reaches 
the saturation temperature somewhat downstream. Vapor 
generated can then exist anywhere in the liquid stream, 
and the bubbles can then begin to coalesce to form the 
slug flow pattern, which progresses into the annular flow 
regime with increasing vapor quality. 

As the quality increases through the saturated 
boiling region, a point is reached where the principal 
mechanism changes from one of boiling, to one of evap- 
oration at the liquid-vapor interface. This is due 
primarily to the establishment of a stable annular climbing 
liquid film and a vapor core, with or without entrained 
liquid droplets. The effective thermal conductivity of 
the thin liquid film on the tube wall is so high that 
insufficient wall superheat exists to allow further 


nucleation. Since boiling is suppressed, this region 
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FIGURE 2 Heat Transfer Conditions for Forced Convection 
Boiling in a Tube 
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is called the two-phase forced convection region of heat 


transfer. As the vaporization proceeds, the heat transfer 
coefficient increases with quality, and very high values 
may be produced, with Nusselt numbers of order 1000 common. 

The transition from principally boiling to principally 
evaporating is complicated and is highly dependent on 
both the flow and heat flux parameters. A further dis- 
cussion will follow in Chapter 2. 

At some critical value of quality the complete 
evaporation of the liquid film occurs, resulting in a 
sharp rise in the wall temperature. The wall is then 
only intermittently wetted by entrained liquid droplets. 
Past the dryout point, the vapor may become considerably 
superheated. 

There has been relatively little work specifically 
in the area of forced convection boiling inside a duct, 
with net vapor generation. Since this is the primary 
region of interest, the next section will present a 
review of some of the better known work. 

1.2 Previous Work in Forced Convection Boiling of Saturated 
Water 

Mumm (1) measured local heat transfer coefficients 
at the exit of a 0.465 in. inside diameter electrically 
heated stainless steel tube for net boiling of water with 
exit qualities up to 70%. The data for qualities less than 


40% were correlated by the equation: 
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> 1.64 
Nu = m + 0.00054 = x|9/^ p90 808 (1.1) 
Vi Arg 
Dengler (2) in 1952 obtained local heat transfer 
coefficients for water in vertical upflow through a i inch 
inside diameter tube, heated by steam jackets. Heat 
transfer coefficients were measured for exit qualities 
up to 80%. It was postulated that the heat transfer 
in forced convection boiling is influenced by the usual 
forced convection effect and a nucleate boiling effect 
when it occurs. At high qualities an increased forced 
convection effect due to the high vapor velocity suppresses 
the nucleation. The rapid increase in volumetric fraction 
of the vapor for small increases in mass quality at low 
and moderate pressures stabilizes the annular flow regime 
even at low qualities. Dengler proposed the following 


correlation: 


h e -0,5 
n = 3.5 (X,,) Fn » (1.2) 
nio tt D 


wher hio is the nonboiling heat transfer coefficient for 
the liquid at the local state and the same total mass 
flow rate obtained from the Dittus-Boelter equation. 


nio * 0.023(51) (08)9- (pr, )9^ (1.3) 
πα 


Xit is the Martinelli parameter widely used for the 


correlation of pressure drop datas 





ai 


Kee (t: Ta) G=)” -— 


Fo represents the correlation factor for conditions 
where nucleate boiling exists. 

In 1960, Sani (3) presented data for water in vertical 
downflow in an electrically heated 0.7194 inch inside 
diameter tube. He compared his data with existing cor- 
relations, but found the deviations substantial. 

In 1962, Schrock and Grossman (4,5) published a 
report on forced convection vertical upflow boiling of 
water in several electrically heated tubes. They derived 
an expression similar to Dengler's for the two-phase 


forced convection region, but introduced the boiling 





number, Bo’ as the extra variable to account for the 
heat transfer enhancement of nucleation. Their correlation 


is of the forms 


n - -0. 67 
3 7400 B, * 0.00015(X,,) (1.5) 


where X.+ and hio are defined as before, and 
° Gh 
fg 
Wright (6) followed Sani's work in downflow boiling, 
and extended the range of flow and heat flux parameters. 
Bertoletti et al (7) obtained heat transfer data 
for steady and transient conditions for vertical upflow 


of water at 1000 psia in electrically heated tubes. Although 








22 
their experiments attempted to determine the critical 
heat flux for a variety of conditions, the behavior for 
heat fluxes less than critical was also explored. 

Chen (8) recognized that there was little consistency 
among the correlations that he examined for water and 
organic fluids, and that none were satisfactory for 
general use. He then proposed a new correlation which 
proved very successful in predicting the forced convection 
boiling heat transfer data he examined. In particular, 
for the water data of Dengler, Schrock and Grossman, and 
Sani, he was able to significantly improve upon the 
correlations proposed by the investigators for their 
own data. This is summarized in Table 1. 

Chen's correlation covers both the saturated nucleate 
boiling region and the two-phase forced convection region. 


It was assumed that both mechanisms occur to some degree 





over the entire range, and that the contributions were 
additive. This method of superposition was first proposed 
by Rohsenow (9). Chen argued that the total heat transfer 
coefficient could be represented by: 


h = haa (1.6) 


ο i hnic 
It was assumed that the convective component, ose could 
be correlated by the Dittus-Boelter equation: 


> 0.85. 0.4 κ 
hac * 0.023Reqp "Prop IP, (1.7) 
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where the Reynolds number, Prandtl number, and thermal 


conductivity are associated with the two-phase fluid. 
Since the heat is ultimately carried through a liquid film 
in annular flow, Chen used the liquid property values. 


Chen defined a factor, F, such that: 
0.8 0.8 





Re 
— = Renp (1.8) 
1 GD(1-x) 
A] 


Equation (1.7) then becomes: 


k 
0.8 0.# ol 
n = 0.023Re, Pry D E (1.9) 


The only unknown factor is the expression for F, a flow 
parameter only, and which Chen suggested would be a 
function of the Martinelli parameter. 

Chen modified the pool boiling analysis of Forster 
and Zuber for the evaluation of hie? the nucleate boiling 


component. Chen proposed the following for the boiling 


contribution: 
y, 0-79, 0.45, 0.59. 0.25 
haig = 0.00122 Íl ] ] So (1.10) 
0.5, 0.29, 0.24 0.24 
S A fg 0v 
0.24 0.75 
AT, AP, 


Because both in pool boiling and forced convection boiling 
the actual superheat is not constant, but decreases with 
distance from the wall, the mean superheat, AT, , in which 
the bubble grows, is less than the wall superheat, AT -- 
For pool boiling, this difference is small, but Chen 
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postulated it could not be neglected in forced convection 


boiling. He then defined a suppression factor, S, such 
thats 


0.99 
s = Δι, (1.11) 


_ Tsat 


Using the Clausius-Clapeyron relation, Equation (1.11) 
can be rewritten: 


.2 0.75 
S i) (se ) ! (1.12) 
Tsat A Psat 


and Equation (1.10) may then be rewritten ass 


10.79, 045, 0.49, 0.25 
haie * 0.00122 ἃς Cpl (1 £o 
κ; 29. 0.24. 0.2 





0.24 = (1.13) 
^Tsat AP at Š 
sat (# T? sat 


Chen further postulated that S could be represented 
as a function of the two-phase Reynolds number, and would 
approach unity at low flow rates and zero at high flow 
rates. He then determined the F and S functions empirically 
from experimental data using an iterative procedure to 
obtain the best fit. These functions are shown in Figures 
3 and +. 

Chen's correlation is at present the best available 
for the saturated forced convection boiling region (10). 
While it is based on physical reasoning with respect to 
the way in which the forced convection and boiling 
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K | Average % Deviation For Correlations 
Data Dengler Schrock & Chen 
Grossman 


Dengler 30.5 20.3 14.7 
89.5 20.0 15.1 















Schrock & 
Grossman 












26.9 48.6 8.5 


Sani 





TABLE 1 


Comparison of Correlations with 
Chen's Data Base 





Reynolds Number Factor, F 
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contributions are superposed, the F and S functions are 


purely empirical relations. In the case of forced con- 
vection, a number of analytical attempts have been made 

to describe the fluid mechanics and resulting heat transfer 
that occurs with an annular flow model (11,12). Additionally, 
considerable work in theoretical models has been done in 
condensation research (13,14,15,16), and although the 

phase distributions may be considerably different for 
condensation, the theoretical models should give results 


that are easily corrected to fit any particular data set. 


1.3 Purpose and Basic Assumptions of Analysis 


It is the purpose of this study to develop a cor- 
relation for the saturated nucleate boiling regions 
using a different and less empirical approach than Chen 
for superposing the effect of the two components. A 
comparison will also be made with Chen's original water 
data, as well as additional data sets not specifically 
used by him in the development of his correlation. 
In particular, the region of interest will be defined 
by: 
1. Saturated, two-phase convective flow of water. 
2. Vertical, axiíial flow in round tubes. 
3. Stable flow. 
h. No slug flow. 
5. No liquid deficiency. 
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6. Heat flux less than critical flux. 


Typically these conditions correspond to annular flow 
or annular flow with entrainment at low and moderate 


pressures, in the quality rage of 1% to 70% (8). 
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Chapter 2 


DEVELOPMENT OF THE CORRELATION 


2.1 General Approach to the Analysis 


In the previous chapter, the correlation of Chen, 
involving a forced convection and a nucleate boiling 
component was shown to be very successful in predicting 
two-phase heat transfer data. The overall approach used 
in this analysis will be similar to Chen's, i.e., an 
additive superposition technique to account separately 
for the forced convection and boiling, but where the 
effect of convection on the point of incipient boiling 
is accounted for. 

The method requires the separate consideration of 
four aspects of the total heat transfer mechanism: 

1. Selection of the forced convection correlation. 
2. Location of the incipient boiling point. 


3. Selection of the nucleate boiling curve 
to be used. 


4. Proper addition of the convection and 
boiling contributions. 
2 Forced Convection Contribution 
Traviss, Rohsenow, and Baron (17) developed a cor- 
relation for forced convection condensation inside tubes. 
Àn annular flow model was used, with the momentum-heat 


transfer analogy using the von Karman universal velocity 





DI 


distribution to describe the liquid film. 


core was assumed to be very turbulent, and the temperature 


in the vapor core and at the liquid-vapor interface was 


assumed equal to the local saturation temperature. 
order of magnitude analysis and non-dimensionalization 
of the heat transfer equations resulted in a simple 


formulation for the local heat transfer coefficient. 


final correlation is of the form: 
0.9 
hpg = F(X,.)Re; Pr, kı 
F^ D 
where Re, = GD(1-x) 
1 


F(X,,) = 0.15(1/X,, + 2.85/X44 


and X.+ is the Martinelli parameter. 


F, = 5Pr, + 51n(1+5Pr, ) + 2. 51n(0.00313Re, 


Their 


(2.1) 


(2.2) 


(2.3) 


0.812, 


(2.4) 


The correlation was very successful in predicting 


condensation data for Freon 12 and 22. 


It was felt that 


such à correlation based on the theoretical analysis of 


the annular flow model would be useful for evaporation 


data as well, although some correction might be required. 


Collier (10) noted that the heat transfer coefficients 
predicted from other analyses for evaporation similar 


to Traviss' were higher than actually observed. Hewitt 


and Hall-Taylor (18) recommended that the heat transfer 


coefficient used for design be about 30% less than that 
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predicted from the theoretical models. 


In order to determine how much, if any, the Traviss 
forced convection correlation overpredicts actual data, 
it is necessary to compare it with data for which there 
is no nucleate boiling. In order to accomplish this, 
the development of an incipient boiling criteria is 


required. 


2.3 Incipient Boiling Criteria 


It is well accepted that in boiling systems, bubbles 
originate at nucleation sites corresponding to cavities 
on the heated surface in which vapor or other gases have 
been trapped. If these sites are modeled as simple conical 
cavities, a bubble will pass through a hemispherical state 
with a radius equal to the cavity mouth radius as it 
grows. For this bubble of vapor to exist, three equilibrium 


conditions must be met at the interfaces 


1. Mechanical (P, = P,) = 27 (2. 5a) 
r 

2. Thermal T,= T4 (2. 5b) 

3. Thermodynamic (H-TS),, - (H-TS)4 (2. 5c) 


Since P. is greater than Pi and T is equal to 
Ti, the liquid must be superheated. The minimum vapor 
temperature necessary for the existence of a hemispherical 
bubble is determined by calculating the difference 
T,-T,a+ at the vapor pressure (See Figure 5). This 


vapor superheat may be related along the saturation line 
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with the use of the Clausius-Clapeyron equation which 


relates the temperature and pressure of two phases in 
equilibrium with the latent heat and volume change for 


a change of phase: 


Neg E (2.6) 
Vfg dl sat 


To integrate Equation (2.6), the following assumptions 


will be made: 


1. Vega V; 
2. We R.T 
P 


ES T, 
Then E m Neg at 
: P- Ry ne 
1 sat 
P h 
σε]. EA = +) (2.7) 
2 v Tsat Ty 
Substituting Equation (2.5a) into Equation (2.7) and 
rearranging: 
TTR 
- = v Sat v ln (1 + 
(Toa) E- aL) (2.8) 


Where the properties are evaluated at the local saturation 
conditions. 

Equation (2.8) represents the locus of stable bubble 
radii in a uniform temperature field. If the bubble is 
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Stable Bubble Radii 
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to grow from its hemispherical state at the cavity mouth, 
there must be net heat transfer to the bubble, and the 
temperature in the liquid layer near the wall must equal 
or exceed T, calculated from Equation (2.8). 

The following procedure to predict intiation of boiling 
was developed by Bergles and Rohsenow (19). The temperature 
gradient in terms of the forced convection heat transfer 


coefficient and wall superheat is: 


κα 3T) 


al) ) (2.9) 


= hpo (Ty Tsat 


The temperature profile near the wall may be approximated 


in a linear form by integrating Equation (2.9): 


T. (y) = 1, "(ahy (2.10) 
J. 


For a given set of flow conditions, both the wall 
temperature and temperature gradient increase with heat 
flux. In Figure 6, a series of lines representing the 
temperature distribution near the wall are shown for 
increasing heat flux. Also shown is the equilibrium 
vapor temperature expressed by Equation (2.8), with 
the cavity radius plotted as distance from the wall. 

The postulated criterion is that nucleation is 
initiated when the temperature gradient is tangent to 
the curve represented by Equation (2.8), implying that 
the liquid temperature equals or exceeds the critical 
value required to nucleate a cavity of radius rerit 


over the entire bubble surface. A slight increase in 
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heat flux would activate cavities in the size yy <T<Yo + 
Numerical estimates of the heat flux required to 
initiate the critical cavity radius may be made by solving 


at the point of tangency: 


T(y)=T, (2.11) 
ati = Iy (2.12) 
dy dr 


Davis and Anderson (20) obtained an analytical solution 


to quations (2.11) and (2.12) with the following results: 


0.5 
E 1 Bk (2.13) 
YA; y, 


u (a/A); ran: 
Urs - Teer), =B + z crit (2.14) 
crit T 
(q/A),, = ELT -T a4) (2.15) 
ib ip * sat'ib i 
where pu 2V Teat"fg (2.16) 
Neg 
At the onset of nucleate boiling: 
(q/A);, = hre (Tw Tsat) ib (2.17) 
Solving Equations (2.15) and (2.17) simultaneously yields: 
2 
1 
Tu Tgarin = #B hpç “o 
T 


The finish of the heating surface can influence the 


boiling curve; the following arguments are best discussed 
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with reference to Figures 7 and 8. 

Figure 7 demonstrates the comparison of rough and 
smooth surfaces which have low non-boiling heat transfer 
coefficients (e.g. natural convection). Under these 
conditions, the wall superheats are relatively high 
even at low heat flux levels, and the liquid temperature 
gradients have a shallow slope. As a result, the initial 
point of tangency will occur at a very large cavity size. 
If this cavity contained trapped gas or vapor, it could 
nucleate. The smooth surface has no cavities in this size 
range and does not nucleate. A slight increase in heat 
flux to q/A, would allow nucleation of the largest cavity 
on the smooth surface, but many of the cavities on the 
rough surface would have been activated with the result 
that boiling would be well established on the boiling 
curve. A higher heat flux, q/Aə, would activate more 
cavities on the smooth surface, and would cause the boiling 
on the rough surface to become even more intense. Asa 
result, two distinct boiling curves of similar slope 
would emerge, with the smooth surface curve shifted to 
the right of the rough surface curve. 

For forced convection with high non-boiling heat 
transfer coefficients, high heat fluxes produce relatively 
low wall superheats, and steeper liquid temperature profiles. 
This results in the initial tangency to occur at an inter- 


mediate cavity size which may be present on both rough 
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FIGURE 6 Temperature Distribution Around a Bubble 
Nucleation Site on a Heated Surface 
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and smooth surfaces. Then the incipient boiling point 
would occur at the same heat flux for both surfaces. 
Increasing the heat flux would tend to activate more 
cavities on both surfaces, with the result that a single 
boiling curve would emerge. 

The above treatment assumes a wide range of "active" 
cavity sizes (residual trapped vapor or gas present). 

In high velocity forced convection, the maximum active 
cavity size may be considerably smaller than the maximum 
size obtained from a surface inspection only. Then the 
inception of boiling would occur at higher fluxes than 

that predicted from the initial tangency. One possible 
reason for this is the inability of large, shallow cavities 
to sustain a trapped vapor pocket under high velocity flow. 
In this case an estimate of the largest "active" cavity 
size must be made and substituted into Equation (2.13) 

to determine the required heat flux. 

Davis and Anderson (20) found reasonable agreement 
with experimental data for water and benzene when a 
maximum active cavity size of lu (3.3X10 ft.) was used. 
Brown (21) measured cavity size distributions on various 
experimental and commercial surfaces and found reasonably 
dense populations of cavities (greater than one site/cm”) 
occur only for equivalent radii less than 10u (3.3X107?ft. ) 
(10). It is felt that these two values represent a range 
for the size of the largest potentially active cavity. 
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FIGURE 7 Comparison of Boiling Heat Transfer from Rough 
and Smooth Surfaces, Low h (21) 
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Equilibrium Vapor Temperaturc 
Superheat, Eq. (2.8) 
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FIGURE 8 Comparison of Boiling Heat Transfer from Rough 
and Smooth Surfaces, High h (21) 
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2-4 Fully Developed Boiling Curve 
In the previous section, the position of the boiling 


curve for forced convections was postulated to be essen- 
tiallly independent of the surface finish. In forced 
convection subcooled boiling, with a sufficiently high 
heat flux, the wall temperature becomes essentially 
independent of the convection effect (see Figure 9). 
This concept will also be applied to the two-phase forced 
convection region, and is the same equivalent approach 
as that of Chen. Three boiling correlations will be 
examined in this study: 
1. Rohsenow boiling correlation (9) 
CoAT sat "Co 7t) em pet? (2.20) 
Neg Ah ee sp] -f,) 
where Cor is a constant for a particular fluid-surface 
combination. 


2. Mikic pool boiling correlation (22) 


(q/A)p pu - B,(fAr, , )"*1 (2.21) 
vem gu M sat 


1/24, 17/8 19/8 (m-23/8). (m-15/8) 
where ger! = Ky / 1 / ; 1 / x / V n 5/ 
-11 -1 5/8 
4, (0009-11 πας ενω 


sat 
By is a dimensional constant which depends on surface 
properties and gravity, and m is a constant which is a 


function of the cavity size distribution. For this analysis, 
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Figure 9 Typical Curve for Subcooled Boiling 





Aly 
m will be assumed equal to 2, in which case Equation 


(2.21) becomes: 
Em. Y x, V? 17/8, pe a 
Ah eg AR a s e a 8s 1/8 


HAT" (2.21) 


3. Thom fully developed subcooled boiling correl- 
ation (23) 
" 0.5 -P/1260 
τα, = W(q/A), °?e (2.22) 


where W is a constant reported by Thom as 0.072 
For each relation, the constant must be evaluated 


for the case of two-phase forced convection boiling. 


2. Superposition T u 
The final part of the correlation is the method of 
adding the forced convection and nucleate boiling com- 
ponents to generate the complete heat transfer curve. 
Two requirements for the curve aret: 

1. At wall superheats below that required for 
incipient boiling, the total heat transfer coefficient is 
equal to the forced convection heat transfer coefficient. 

2. At sufficiently high superheats, the heat 
transfer coefficient approaches the fully developed boiling 
curve. 

The procedure to be used in this analysis is to force 


the boiling component to be zero at the incipient boiling 
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points 
q/A 7 q/Ag, * q/Ag - q/A' (2.23) 
where q/A' is the heat flux on the fully developed boiling 
curve at the incipient wall superheat calculated from 
Equation (2.15) or (2.19). 
If the fully developed curve has a slope of non 


log-log coordinates, thens 


n 
q/A * q/Agc * q/Ag(1 Cea) Jos (2.24) 


sat 


n 
so that 1 I in) is a suppression factor for nucleate 
T 
sat 
boiling. Figure 10 shows a sample plot of the superposition 


procedure. 
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FIGURE 10 Superposition Technique 
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Chapter 3 


CORRELATION OF DATA 


3.1 Data Base and Properties 


Eight sets of water data taken in round tubes were used 
as the basis for comparison with the correlation. Table 2 
shows the range of experimental conditions covered by the 
data. 621 data points consisting of local conditions of 
mass flow rate, heat flux, saturation temperature, quality, 
and measured heat transfer coefficient were examined. 

The water properties used in the calculations were 
linearly interpolated from the 1963 International Skeleton 


Steam Tables (24) after conversion to engineering units. 


3.2 Forced Convection 


The incipient nucleation criteria outlined in Section 
2.3 was applied to each data point, with the heat transfer 
coefficient required in the analysis calculated from the 
unmodified Traviss forced convection correlation (Equation 
2.1) as a first estimate. Equation (2.18) was used to 
first estimate the incipient heat flux based on the point 
of tangency, from which the size of the critical cavity 
radius was calculated from Equation (2.13). If this cavity 
was larger than the assumed maximum active cavity radius, 


the required heat flux to boil was calculated from: 
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Bk; he 


2 
(q/A),, » “max _ (3.1) 
k; _ h a 
— F 


Tmax 


which is obtained by solving Equations (2.14) and (2.17) 


simultaneously with net a 


Several maximum cavity radii were tried; Pao fie 
was found to give the best results, and was within the lmits 
previously described in Section 2.5. Figure 11 graphically 
displays this resulting nucleation criteria for water. The 
incipient heat flux may also be found by plotting the 
forced convection heat transfer component and reading off 
the value of heat flux or wall superheat at the intersection 
with the nucleation line of the correct pressure. 

Those points with measured heat flux less than the 
incipient heat flux were classified as non-boiling, and 
plotted against the Traviss parameters: 


NuF, 


0.9 
1 


vs. Xit 


Re Pr 


Í 


with Nu = haataP 


Ey 
A curve of the same general form as the original Traviss 


F(X,,) was then fit through the data; i.e. 
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FIGURE 11 Incipient Nucleation Criteria for Water 
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PO py = KI (1/X., + KS 7X7) (3.2) 
The procedure was repeated through three iterations with 
the new expressions for P(X try inserted in Equation (2.1). 
It should be noted that the use of the proposed nucleation 
criteria at this stage is only as an indicator for possible 
boiling, and has no effect on the final equation for F(X,,)-. 
The proposed correlation for the forced convection 


component is: 


0.9 

has = Rei PrjF(X,,) Ky (2.1) 
F, D 
F(X,,) = 0.15(41/X,, + 2.0(1/X,,)°°7*) (3.4) 
F, = 5Pr, + 51n(1+5Pr)) + 2. 51n(0.00313Re, 9:912) 
2.4) 

Re, = GD(1-x) 

M1 


This correlation is compared with the original Traviss 
correlation in Figure 12, and the fit of the non-boiling 
data to the proposed correlation is shown in Figure 13. 

The same non-boiling points were plotted against the Chen 


macroconvective parameters, 


Nu vs. 1/X 


0.4 


0.8 


and are shown in Figure 14. The two correlations were 
compared by computing the average deviation between the 


predicted and experimental values of the forced convection 
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parameters. 


For the Hall-Traviss correlation, 
(F(X = (FCG) 


Deviation = UP red data 
(Ft Naar 


For the Chen correlation, 
- F 


Deviation = ? pred data 
F data 

The results of this comparison are summarized in 
Table 3, and shows that the proposed forced convection 
correlation is generally superior to the Chen macrocon- 
vective correlation, when applied to data for which there 
is no nucleation. It should be noted that the Chen method 
requires the superposition of both a forced convection and 
nucleate boiling component, and that for the typical range 


of interest, the nucleate boiling component is never com- 


pletely suppressed (see Figure 4). 


3.3 Nucleate Boiling 


For those data points where the applied heat flux was 
higher than the incipient flux, enhanced heat transfer due 
to nucleate boiling was assumed to take place. This effect 


can be determined by rearranging Equation (2.24): 


a/a, = (YAgata l q/Apo) (3.5) 
1 - AT sat E 
5 Tsat,dat 
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V/A sata 
data 


and A Tsat,data " 


Knowing the residual heat transferred by nucleation and 

the corresponding wall superheat, the appropriate value 

of C p» By, or W (Equations (2.20), (2.21), (2.22)) can be 
calculated for a particular data point. Following the 
argument for the effect of surface finish on forced con- 
vection boiling, a single point should establish the pos- 
ition of the boiling curve. In reality, the large scatter 
present in the prediction of incipient boiling and nucleate 
boiling phenomena led to a graphical approach to obtain 

a "best" value of the particular constant. In this respect, 
a value judgment was made by weighting those data points 
where significant boiling should exist (low G, low x, high 
q/A). The best values obtained for the boiling constants 


were: 


Cor = 0.0288 


B, = 0.0000213 


M 





F(XTT) = (NU=*FƏ) Z(PRL*REYL*<=0.9) 
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Traviss 0.476 
F(X, 4)=0.15(1/X,,+2.85(1/X,,) ) 


Modified Traviss 0.32 
F(X, ,)=0.15(1/X,,+2.0(1/%,,) 





1055 101 109 10° 
MARTINELLI PARAMETER, XTT 


FIGURE 12 Comparison of Traviss (17) Forced Convection 
Correlation and Proposed Forced Convection 
Correlation 





F(XTT) +» (NU*F2)/(PRLAREYL+*O, 9) 
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P(X,¿)50.15(1/X,,+2.0(1/x,,)9:32) 


< 
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FIGURE 13 Comparison of Proposed Forced Convection 
Correlation with Non-Boiling Data 





F += NU/ (O, O23*REYL +0, B»PRL +0, 4) 
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FIGURE 14 Comparison of Chen (8) Forced Convection 
Correlation with Non-Boiling Data 
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Data # 


Dengler 25 


Schrock & Grossman 9 
Series A 


Schrock & Grossman O 
Series E 


Schrock & Grossman O 


Series F 

Bertoletti 1 
Sani 83 
Wright Series 1 67 
Wright Series 2 29 
Average 214 


A 
072677 


0.1310 


0.5514 
0.0886 
0.0895 
0.1867 
0.1270 


Forced Convection Correlations 


A - Modified Traviss Forced Convection (Equations 


ο. ος ο, 9.2) 


Average Deviations for 
Non-Boiling 


B 
0.3626 


0.1531 


0. 5238 
0,1252 
0.1188 
021772 
0.1610 


B - Chen Macroconvective (Equation 1.9) 


TABLE 3 


Comparison of Forced Convection Correlations with 


Non-Boiling Data Predicted by Proposed Incipient 
Boiling Criteria 
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Chapter 4 
RESULTS OF THE ANALYSIS 


With all four components of the correlation complete, 
it was tested against the data. This method may assume 
either the wall superheat or the heat flux specified. Since 
a specified heat flux is the more usual design case, the 
data was reduced using this assumption, entailing an it- 
erative procedure to determine the wall temperature re- 
sulting from the experimental heat flux. Sample calculations 
for both approaches are provided in Appendix II. Due to 
the large number of data points examined, a digital computer 
was employed, but the predicted heat transfer coefficients 
are readily hand calculated. A listing and sample output 
from the data reduction program are provided in Appendix IV. 

The criteria chosen for comparison between predicted 
and experimental results is the deviation, or difference 


ratio. 


For the case of specified wall superheat: 


Deviation - (n red I haata! 
haata 


For the case of specified heat flux: 


Deviation - PTT ` AT at data’ 


ATsat ‚data 





p 


The results were produced in two forms. The first 
is scatter plots of predicted heat transfer coefficient 
vs. the experimental value (see Appendix III). Points 
lying on the 45 degree line represent perfect correlation. 
The scatter plots show that for the most part, the data 
is roughly centered about the 45 degree line. The anomo- 
lous points well off the experimental values can be found 
on each plot, and are possibly bad data. 

The second presentation of results is a tabular summary 
of the computed average absolute value deviation (see Table 
4), The proposed method of correlation compares very 
favorably with the results from the widely accepted Chen 
method, and in fact, the Hall-Traviss forced convection/ 
Mikic nucleate boiling had the lowest combined average 
deviation of +15.4%. Each of the three proposed correl- 
ations, us ing the methods outlined in Chapters 2 and 3, 
was able to predict heat transfer coefficients over a 
wide range of conditions with an average deviation less 


than +30% for any particular data set. 
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Average Deviations for 
Correlations 
Data Set # Points I IET III IV 


Dengler 119 {0.2036 } 0.2495 | 0.1619 [10.1595 
0.2096 | 0.2034 | 0.2302 | 0.1902 


















Schrock & Grossman 160 
Series A 











0.1886 | 0.1637 





Schrock & Grossman 50 0.1692 | 0.1704 


Series E 






















Schrock ἃ Grossman 38 0.2179 | 0.1442 | 0.2856 | 0.1546 


Series F 








0.1842 





0.2013 
0.0875 
0.0895 


0.1845 
0.0879 
0.0895 


0.2009 
0.0947 
0.0938 





Bertoletti et al 64 









0.0876 
0.0895 





Sani 84 
Wright #1 67 














Wrisht #2 39 0.1638 | 0.1772 | 0.1827 | 0.1812 
Correlations 

I Chen Correlation 

TI Hall-Traviss Forced Convection with Rohsenow 


Nucleate Boiling, C „p70. 0288, r nax” O. 00001 fb 


ΠΤΙ Hall-Traviss Forced Convection with Thom 
Nucleate Boiling, W=0.132, Tax 0.00001 το. 


ΙΥ Hall-Traviss Forced Convection with Mikic 
Nucleate Boiling, By70. 0000213, Pax 0 00001 ft. 


TABLE 2 


S ar Results 
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Chapter 5 
SUMMARY AND CONCLUSIONS 


The case of forced convection boiling of saturated 
water was treated in this study with the hope of developing 
a method of correlating local heat transfer coefficients 
based on physical principals. The total heat transfer 
mechanism was postulated to be made up of a forced convection 
effect and a nucleate boiling effect where it exists. 

For annular flow, which exists over much of the quality 
range for low and moderate pressure, the forced convection 
component is due to the convective transport through a 
highly turbulent liquid film. A modified Traviss (17) 
analysis based on the momentum-heat transfer analogy, 
using the universal velocity profile to describe the annular 
liquid film was used. The simplified design equations 
proposed by Traviss were corrected to fit non-boiling data. 
The incipient nucleation criteria proposed by Bergles and 
Rohsenow (19) based on vapor bubble equilibrium in a linear 
temperature profile near the wall was employed to account 
for the effect of forced convection on the onset of boiling. 
A maximum active cavity radius of 107 ft. was suggested. 
Boiling was assumed to exist only if the applied heat flux 


was greater than the incipient boiling heat flux. 
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The proposed forced convection correlation was found 
to fit the non-boiling data identified by the incipient 
nucleation criteria better than the Chen (8) macroconvective 
relation. 

The nucleate boiling and forced convection components 
were superposed by forcing the boiling heat flux to be 
zero at the point of incipient nucleation. The position 
of the fully developed boiling curve was argued to be 
essentially independent of surface finish using the same 
reasoning as Brown (21). The nucleate boiling correlations 
of Rohsenow, Mikic, and Thom were each tried as the rep- 
resentative component to be added to the modified Traviss 
forced convection component. 

The proposed method was used to predict heat transfer 
coefficients for eight sets of water data, assuming the 
heat flux as the specified independent variable. The 
widely recommended Chen correlation was also applied to the 
same data. The proposed method produced results which 
compared very favorably with the local values predicted by 
the Chen method. The Hall-Traviss forced convection/Mikic 
nucleate boiling had the lowest combined average deviation 
of all correlations tested (+15.4%), while the Chen method 
produced a combined average deviation of +17.4%. The 
results were found to be relatively insensitive to the 


particular boiling correlation selected. 
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The proposed method has the advantage of not requiring 
the use of purely empirical and graphically presented 
correlation factors, and is recommended for predicting 
forced convection boiling heat transfer coefficients for 
water at less than 1000 psia, in the approximate quality 
range of 1-70%. 

Although the same mechanisms of heat transfer should 
apply to other fluids, the coefficients in the equations 


for forced convection and boiling should be verified before 


general use. 
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Appendix I 


PROPERTIES OF SATURATED WATER 


KEY 


C SATURATION TEMPERATURE; ΠῈΣ ΕΛΗΞΕ 

ER SATURATION PRESSURE, PSIA 

ERE TOUID DENSITY, LBM/FT**3 

Bev NrOR DENSITY, LBM/ET**3 

S-LATENT HEAT OF VAPORIZATION, ETU/LPM 
O ID SPESLETC REAT; BTU/L3M-DEGF 
01 SURFACE TENSION, LBF/FT 
ELLIOUID VISCOSITY, EBM/HPR-FT 

ο ΛΡΟΕ VISCOSITY, LBM/*3-rT 

BNS-LIOUID THERMAL CONDUCTIVITY, BTU/HR-FT-DEGF 
11-LIQUID PRANDTL NUMBER 
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APPENDIX II 


Sample Calculations 


Sample 1 - S ) is specified 


sat 


Given: 


D 1.0 in 


2900 lbm/hr 


m 
" O 
AT πο 17.9 “F 


s 


2 


h = 5530 BTU/hr-ft 


data 
x = 0.074 


255 Op 


r3 
N 


sat 


rg 
"n 


sat 33 psia 


Properties: 
(1 
(v 


he 


58.4 1bm/ft? 
0.080 lbm/ft2 


en 941.7 BTU/lbm 
M, = 0.54 1bm/hr-ft 
4, = 0.031 1bm/hr-ft 
p1 = 1-02 BTU/1bm-°F 
Pr, 1.38 


k, = 0.397 BTU/hr-ft-?F 


ς 


S= 0.0037 1b/ft 
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I. Forced convection heat transfer coefficient 


Hall-Travi 


0.1 0.5 
Calculate X.. acp 1) E 
My 1 


is E -9 0.080 9.3. Yo? 
O. 074 (9: 031 





G=m= 2900 = 531700 1bm/hr-ft* 
A 0.0055 


P(Xp4) » 0.15( 1/X,, * 2.0(1/x,,)0* 22) 
= 0.697 


Re 5Pr, $ 51n(1+5Pr, ) + 2. 51n(0.00313Re ,9* 912) 


2 


Re, = GD(1-x) = 700 1119220) = 76000 
Mi 0. 54)(12 


P, = 5(1.38) + 51n(1+5(1.38)) + 
2. 51n(0.00313(76000)9* 912) 
P, = 25.63 
0.9 
_ Re χει de 
he, = Rey SL a 


F2 


6000) 020.697) (1.38) 0.397(12) 
| 25.63 1 


4426 BTU/hr-ft^-9F 
Chen Forced Convection 


1/X,, = 2.11 
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72 
F= 4.5 (Figure 3) 


"mac 


A 


0.8 0.4 
0.023 Re, Pry TOEF 


D 
E 0.8 0.4 
= 0.023(76000) (1.38) (0.397) (12) (4. 5) 
1 
2 0 

= 4500 BTU/hr-ft”-"F 
Incipient Boiling Point 
Calculate 

Benen sat’ =" C «00 255+460) (12.48) 
Ney 941.7)(778 


9.01X107? ft?R 


= £ => 
AT ot) ap = A Bhac - a. το ) (4426) 
Nee : 


4.0 Ἢ (based on tangency point) 


Bk 
T -a = | 1 


(q/A);+ = heal Toot) ap = (4426)(4.0) = 17700 


i 


Pomit =f (201X107?) (0.397) = 4.5X107> ft. 
17700 
Dort is greater than poe SO D uec. ix 
EU 
P ns 10 ft. 
ER. 
Απο -  Ὁ : = 10.1 °F 
sat'ib ` E. ° 
τ 
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III. Boiling Heat Transfer Coefficients 


Chen 





Re,F? 53 - (76000) (4.5)1*79 & 5.0Xx10° 


S = 0.1 (Figure 4) 


0.79, 0.45, 0.49, 0: 25 
h... = 0.00122 Kı “CP; T 


mic 
0.5. 0.24 0.24 D; 24 
y My Neg ἐν 


0.24 0.75 
x AT at APsat S 
AP „+ = Mp fg AT, an (Clausius-Clapeyron) 
Tsat” fg 


hie = 0.00122(0.397)9* 79(1.02)9*^5(58, 49:9 


E (0.0037)%3(0, 5439-24 (911 .7)9.28 


8.0.25 D 24 o 
X (4.173X10°) AT, AB as 


(0. ogo)9 25 


he Ne, AT sat 
AP 41 8 AT. 
sat pue μας fey 


sat" er 


AP 


es 82.11 AT. at 


0.75 = 0.75 
AP = 27.28 AT, „+ 


sat 


h 


11.36 E = 11.36(17.9)°° 77 


198 BTU/hr-ft^-?F 


mic 


Amic 


Thom Nucleate Boiling will be used to demonstrate 


proposed method of correlation 
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n 
Ur reg = q/Apo + q/Ap(1 - (Aisat.in) ) 


sat 


= 2 for Thom correlation 





P/1260\2 
q/As = “sate 
| 33/1260 
= (12: e 
0.132 
q/A (1 - eM - 
Tsat 
h, = WAR = 13209 = 7 
NP. 17.9 


(W = 0.132) 


2 
| - 19378 BTU/hr- 


19378(1 - (19:05 
17.9 


13209 BTU/hr-ft* 


2 0 


38 BTU/hr-ft^-?r 


Total Heat Transfer Coefficients and Deviations 


5520 


ft? 


Chen 
hs pe: 4r Ami = 4500 + 198 = 4698 BTU/hr-ft 
Deviation = κ. _ Naata = 4698-5530 
Nagata 5530 
= -0.15 
Hall-Traviss/Thom 
h = hp, + hy = 4426 + 738 = 516b BTU/hr-ft^-? 
Deviation = 5164-5530 = -0.067 


2 0 


E 


F 
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Sample 2 - g/A is specified 


Given: 


Identical conditions as Sample 1, with 


2 


exception that q/A = 99000 BTU/hr-ft”, 


data 
and the wall superheat must be predicted. 


I, Forced Convection Heat Transfer Coefficient 


Chen ha. = #500 BTU/hr-ft“-9F 


Hall-Traviss hag = 4426 BTU/hr-ft“-9F 


II. Incipient Boiling Point (not applicable to Chen method) 
q/A;, * 4700 BTU/hr-ft^ 

E ο 

GT atip # 10.1 F 


III. Boiling Heat Transfer Coefficients and Deviations 


Chen 
_ q/A _ ο 
Let A Tsat/try = c = 22 “F 
mac 
! 0.99 . 0.99 
ης... = 11.36(AT. op) try = 11.36(22) 
= 242 
heey = 4500 + 242 = 4742 
(AT +) = Q/Aaata = 99000 = 20-9 °F 
hee 22 
o ο 
Now let io. at ry = 20.9 ?F 
- 0.99 _ 
νυν = 11.36(20.9) = 230 
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Mery = 4500 + 230 = 4730 BTU/hr-ft 


O 
πο UY 5 που (£ 


Assume the trial and error solution has 


2 Op 


converged to the values 


AT = 20.9 °F 


sat 
h = 4730 BTU/hr-ft^-?F 


) - (ÔT a4) 


Deviation = AT d tl data 
TT 


sat'data 


= 20, = = 0.168 
1 


Hall-Traviss Forced Convection/Thom Boiling 


2 
33/1260 AT ες) 2 
q/Az = (122/1250 Asa: ) = 60.47 AT t 


0.132 
Let (AT iu, 7 q/Aqata = 99000 = 22.4 °F 
y hs 4426 


2 
q/A,(1 (BT at. ib ) (60.87) (22.4) (1- (10.11) 
AT 22.4 





sat 
= 24200 BTU/hr-ft? 
ha = 9/Az = 24200 = 1080 
AT ot) try 22.4 


h = hec + hg = 4426 + 1080 = 5506 
_ q/A a my ο 
(OT atl try = DU = 22006 = 18.0 “F 


EN O 
Iu * 18.0 °F 





Tu 


Following the above procedure, the solution 


will converge in four iterations to the values: 


_ ο 
E - 19.9 F 
h = 5243 BTU/hr-ft^-^F 
Deviation = 18.9 - 17.9 = 0.056 


17. 
17.9 
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Appendix III 


GRAPHICAL COMPARISON OF CORRELATIONS WITH DATA 
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Data of Schrock and 
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H, EXPERIMENTAL BTU/(HR-SQFT-DEGF) 





10 





BTU/ (HR-SQFT-DEGF) 


H, PREDICTED 


Data of Wright, Series 2 


Hall-Traviss FC 
Thom Boiling 


10^ 


H, EXPERIMENTAL BTU/ (HR-SQFT-DEGF ) 
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Appendix IV 
DATA REDUCTION PROGRAM AND SAMPLE OUTPUT 


A data reduction computer program was written to 
expedite the analysis of the large number of data points. 
The program requires the input of local conditions of mass 
flow rate, heat flux, vapor quality, saturation temperature, 
and experimental heat transfer coefficient. For each 
of the three proposed as well as the Chen correlation, the 
program uses the heat flux as specified, and predicts the 
resulting wall superheat using an iterative trial and error 
method, superposing the forced convection and boiling 
components. 

For each data point, the program outputs the mass 
velocity, heat flux, saturation pressure, saturation temp- 
erature, quality, experimental wall superheat, experi- 
mental heat transfer coefficient, Martinelli parameter, 
liquid Reynolds number, Hall-Traviss forced convection 
heat transfer coefficient, incipient boiling heat flux, 
and the predicted heat transfer coefficients and resulting 
deviations for the Chen, Hall-Traviss/Rohsenow, Hall-Traviss/ 
Mikic, and Hall-Traviss/Thom correlations. 

Water property values are determined using a linear 


interpolation table search routine. The empirical Chen 
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F and S factors are determined in a similar manner. 


The job control language, input/output device numbers, 


and plotting routine are unique to the MIT Interdata Joint 


Computing Facility, but the remainder of the program is 


standard Fortran IV. 


The organization of a sample input deck is shown below: 


// XEQ 3 
BC/1000/ 
LC/AREA 1/1000 
LC/AREA 2/1000 


CARD 1 
CARDS 2-5 


CARDS 6-9 


CARDS 10-13 
CARDS 14-17 


CARDS 18-21 


CARDS 22-25 
CARDS 26-29 


CARDS 30-33 


X-axis label for scatter plot, col.1-40 
Y-axis label for scatter plot, col.41-80 
FORMAT 40A2 


25 values of saturation temperature (°F), 
FORMAT 8F10.6 


25 corresponding values of liquid density 
(1bm/ft2J), FORMAT 8P10.6 


i y τ A 2 vapor dens ity 


> " T latent heat 
(BTU/1bm), FORMAT 8F10.6 


= n S 2 D liquid specific 
heat (BTU/lbm-°F), FORMAT 8F10.6 


S ü B E u liquid surface 
tension (lbf/ft), FORMAT 8F10.6 


" 1. 5 5 E liquid viscosity, 
(1bm/hr-ft), FORMAT 8F10.6 


" " " " " — vapor viscosity 
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CARDS 34-37 25 corresponding values of liquid thermal 
conductivity (BTU/hr-ft-°F), FORMAT 8F10.6 


CARDS 38-41 E = š N " Prandtl number 


CARDS 42-45 ° n š : "saturation press- 
ure (psia), FORMAT 8F10.6 


CARDS 46-49 28 values of 1/X 44 in the range 0.1 to 100 
FORMAT 8F10.6 


CARDS 50-53 28 corresponding values of Chen F-factor 
FORMAT 8F10.6 


CARDS 54-56 21 values of Chen two-phase Reynolds number 
in the range 10° to 10’, FORMAT 8F10.6 


CARDS 56-58 21 corresponding values of Chen s-factor 
CARD 59 Inside tube dimater (in.), maximum active 


cavity radius (ft.), Cof’ Bpr Mm W, number 
of points in data set FORMAT 6E12.4,18 


CARD 60 Name identifier of data set (col.1-80) 
FORMAT 4OA2 
CARD 61 Fluid identification (col.1-80) FORMAT 40A2 
CARD 62 Flow orientation (col.1-80) FORMAT 40A2 
CARD 63- Mass flow rate (lbm/hr), heat flux (BTU/hr-ft^) 


62+#pts. , quality (decimal fraction), saturation temp- 
erature (°F), experimental heat transfer co- 
efficient (BTU/hr-ft^-?F) FORMAT 5F15.4 
CARD 62+#pts. Blank card if this is only data set or re- 
peat with card 59 for batch runs of data 


sets; a blank card must follow the last 
data point of the last data set 


The program listing and output from the calculations 


are provided in the following pages. 
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